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Abstract 

We study the large-time behavior of solutions to the compressible Navier-Stokes equations 
for a viscous and heat-conducting ideal polytropic gas in the one-dimensional half-space. A 
rarefaction wave and its superposition with a non-degenerate stationary solution are shown to 
be asymptotically stable for the outflow problem with large initial perturbation and general 
adiabatic exponent. 
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1 Introduction 

The one-dimensional motion of a compressible viscous and heat-conducting gas in the half space 
]R_i_ := (0, oo) can be formulated by the compressible Navier-Stokes equations 

pt -I- {pu)x = 0, 

< {pu)t + {pu^ + P)x = {pux)x, ( 1 - 1 ) 

^ {pE)t + {puE + uP)x = (kOx + PUUx)x, 

where t > 0 and x G R+ stand for the time variable and the spatial variable, respectively, and 
the primary dependent variables are the density p, the velocity u and the temperature 0. The 
specific total energy E = e -I- with e being the specific internal energy. It is known from 
thermodynamics that only two of the thermodynamic variables p, 0, P (pressure), e and s (specific 
entropy) are independent. We focus on the ideal polytropic gas, which is expressed in normalized 
units by the following constitutive relations 

P = Rp9, e = c„0, s = Cy ln(p^~'’'0), (1.2) 
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where i? > 0 is the gas constant, 7 > 1 the adiabatic exponent and c„ = i ?/(7 — 1 ) the specihc 
heat at constant volume. Positive constants fj, and n are the viscosity and the heat conductivity, 
respectively. 

The system (1.1)-(1.2) is supplemented with the initial condition 

(p, u, 0)\t=o = (po, uo, 6 »o), (1.3) 

which is assumed to satisfy the far-held condition 


lim (po, Uo, do){x) = {p+, u+, e+), (1.4) 

x—^oo 

where p+ > 0, u+ and 0-|_ > 0 are constants. For boundary conditions, we take 

(■u,6>)(t,0) = (u_, 6 »_), (1.5) 


where u_ and 6- > 0 are constants. The initial data (1.3) is assumed to satisfy certain compati¬ 
bility conditions as usual. 

The boundary condition u(t, 0) = u_ < 0 means that the huid blows out from the boundary, and 
hence the initial boundary value problem (1.1)-(1.5) with u- < 0 is called the outhow problem. The 
problem (1.1)-(1.5) with u- = 0 is called the impermeable wall problem, which has been studied 
in [ 6 , 7, 20, 21, 31] and so on. According to the theory of well-posedness for initial boundary value 
problem, one has to impose one extra boundary condition p(t, 0 ) = p_ on {x = 0 } for the case 
when u- > 0. This case is called the inflow problem and has been investigated by Matsumura 
et al. [4, 6 , 9, 22, 27, 28]. We refer to Matsumura [19] for a complete classihcation about the 
large-time behaviors of solutions to initial boundary value problems of the isentropic compressible 
Navier-Stokes equations in the half space R+. 

The main purpose of this article is to study the large-time behavior of solutions to the out¬ 
flow problem (1.1)-(1.5). The nonlinear stability of the stationary solution, the rarefaction wave 
and their composition has been addressed in [15, 26] under small initial perturbation. For large 
perturbation case, Qin [26] proved that the non-degenerate stationary solution is asymptotically 
stable under the technical assumption that the adiabatic exponent 7 is close to 1. Recently, Wan 
et al. [30] establish the asymptotic stability of the non-degenerate stationary solution with large 
initial perturbation and general adiabatic exponent 7 . In this article we are going to study the 
case when the corresponding time-asymptotic state is a rarefaction wave or its superposition with 
a non-degenerate stationary solution under large initial perturbation. 

We first investigate the large-time behavior of solutions toward the rarefaction wave for the 
outflow problem or the impermeable wall problem (1.1)-(1.5). To this end, we assume that positive 
constants p+, u± and 6± satisfy 

^{e_/fl+)Tip+ . - - - 

u-=u++ Y ^7P+ 2 dz < u+, (1.6) 

Jp+ 

so that (p+,it+, 0 +) G Rz{p-,U-,9-) for 

P- = P+- (1-7) 


Here R^[p-,U-,6J) is the 3-rarefaction wave curve through {p-,u-,0-) given by 
i? 3 (p_,u_, 6 »_) := \ {p,u,0) 


p > P-, p^ ^0 = p[_ ^0-, 


u = U--\- y RlP- ^0-z 2 dz I 

J p. 
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We assume further that 

u_ + y'R'yB- > 0. (1.8) 

Then as time t goes to infinity, the solution of the problem (1.1)-(1.5) is expected to converge to the 
3-rarefaction wave {p^,u^,d^){t,x) connecting (p_,u_,0_) and (p+,tt+,0+), which is the unique 
entropy solution of the Riemann problem for the corresponding hyperbolic system of (1.1)-(1.2) 
(i.e. the compressible Euler system) 

{ Pt + {pu)x = 0, 

{pu)t {pu^ -f P)x = 0, (1-9) 

{pE)t + {puE uP)x = 0 


for (t, x) € M+ X K with initial data 


(p,u,0)(O,x) 


(p-,U-,0-) for a; < 0, 

(p+,M+,0-i_) for a; > 0. 


( 1 . 10 ) 


Note that {p^,u^,9^){t,x) = {p-,u-,6-) for each {t,x) G [0,oo) x (—oo,0] due to (1.8). 

We construct a smooth approximation {p,u,6) of {p^,u^,9^) for deriving the stability of the 
rarefaction wave. As in [9], we consider the Cauchy problem for the Burgers equation 


Wt -f wwx = 0 for (t, x) G IR+ X 

f-X + 

w{0,x) = W- -I- wkg / z'^e~^dz for x G K, 

Jo 


( 1 . 11 ) 


where w := w+ — w-, g > 16 is some fixed constant, x+ := max{x, 0} is the positive part of x, 


and the constant kq satisfies 


kq / z'^e ^dz = 1. 


The smoothed rarefaction wave {p,u,9) connecting {p-,u-,9-) and (p+,u+,9+) is defined by 
X3{p,u,9)(t,x) = w{l -\-t,x), 

{p^-^9){t,x) = 

u{t^x)=u+-\- / y R'ypj^ ^9^z = dz, 

J pi 


( 1 . 12 ) 


where p- is given by (1.7) and w(t,x) is the unique solution of (1.11) with w± = X3{p±,u±,9±). 

Now we state our stability result of the rarefaction wave (p^, u^, x) to the outflow prob¬ 

lem or the impermeable wall problem with large initial perturbation. 

Theorem 1. Assume that (p+,u±,0±) and the initial data (pq,uo,6>o) satisfy (1.6), (1.8) and 

inf {po(a;),6>o(x)} > 0, {po - p,uo - u,9o - 9) e H^{R+). (1.13) 


Then there exists a positive constant ei such that if u- < 0 and the boundary strength S := 
|(tt+ — U-,9+ — 0_)| < ei, the initial boundary value problem (1.1)-(1.5) admits a unique solution 
(p, u, 9) satisfying 

{p — p,u — u,9 — 9) G ^([0,oo);i7^(K+)), 

Px- Px & T^(0,oo;L^(M+)), {ux - Ux,9x - 9x) G L^(0, oo; i7^(IR+)). 


(1.14) 
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Furthermore, the solution {p,u,0) converges to the rarefaction wave {p^, 
tends to infinity: 



, 9 — 9^){t, a;)| = 0. 


u^, 9^) uniformly as time 


(1.15) 


Remark 1.1. For the Cauchy problem to the compressible Navier-Stokes equations (1.1)-(1.2) with 
generic adiabatic exponent 7 in the whole space R, we can employ the methodology developed 
in this paper to obtain the time-asymptotic stability of the rarefaction waves under large initial 
perturbation, which extends the corresponding stability results in [14, 25] for the case with small 
initial perturbation or the case when 7 is close to 1. We refer to a recent work [11] for the stability 
of superposition of viscous contact wave and rarefaction waves to the Cauchy problem for the 
compressible Navier-Stokes equations in Lagrangian coordinate. 

Next we intend to study the time-asymptotic stability of the superposition of a non-degenerate 
stationary solution and a 3-rarefaction wave. For this purpose, we let (p, u, 9) be the stationary 
solution of (1.1)-(1.5) connecting (u-,9-) and (pm, Um, 9m), namely (p, u, 9) depends solely on the 
variable x and satisfies 

( {pu)' = 0, 

J {piif + Py = pu", (1.16) 

[ {puE uPy = k 9" + p^uu'y 

for X S M+ and 

(u, 0)(O) = (m_,0_), \im {p,u,9){x) = {pm,Um, 9m)- (1-17) 

X—¥00 

where P := Rp9 and E := c„0 -I- . A stationary solution (p, u, 9) is called to be non-degenerate 

if for each n G N, 

\d)f{p- Pm,u- Um,d - 9m){x)\ < CSe~°^, (1.18) 

where C, c are positive constants and S := \{um — U-,9m — 9-)\ is the boundary strength of the 
stationary solution. The existence of (non-degenerate) stationary solutions has been shown by 
Kawashima et al. [15] and will be restated in section 2. 

We assume that 


(p-|_, ri-j-, 9^) G R^ i^Pm, a-m, 9m), \lR'^9m ^ Um ^ A9) 

so that there exist a 3-rarefaction wave {p^,u^,9^) connecting {pm,Um,9m) and (p+,M+,0+). It 
is expected that the large-time behavior of solutions to the outflow problem (1.1)-(1.5) is deter¬ 
mined by the composition (p, u, 9) of the stationary solution (p, u, 9) and the 3-rarefaction wave 
(p^, u^, 9^): 

(p, U, 9){t, x) = (p, U, 9){x) + (p^, U^, 9^){t, x) - (pm, Um,9m)- (1-20) 

In order to derive the stability result, we introduce the smoothed asymptotic state 

(p, U, 9)(t, x) = (p, U, 9)(x) + (p, U, 9)(t, x) - (pm, Um, 9m), (1-21) 

where (p, u, 9) is the smooth rarefaction wave connecting {pm, Um, 9m) and (p+, 0_|.). We define 

S := \{um - u+,0m - 6 *+)]. 

Under the above preparation, we have the following stability result on the superposition (p, u, 9) 
with large initial perturbation. 
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Theorem 2. Assume that there exists a non-degenerate stationary solution {p, u, 6) connecting 
{u_,9_) and Assume further that (p_|_,u+,0+) and the initial data {po,Uo,9q) satisfy 

(1.19) and 

inf {po(a;), 6 >o(a::)} > 0, {po - p,uo - u,9o - 9) G (1.22) 

Then a positive constant e 2 exists such that the outflow problem (1.1)-(1.5) with S -\- 5 < €2 has a 
unique solution {p, u, 9) satisfying 

{p- P,u-u,9-9) G C{[0,oo);H^{R+)), 

Px- Px& L^(0, oo;L^(M+)), {u^ - Ux,9^ - 4) G L^(0, 00 ; i/^(IR+)). 

Furthermore, the solution [p, u, 9) converges to the composition (p, u, 9) of the stationary solution 
{p,u,9) and the rarefaction wave {p^,u^,9^) uniformly as time tends to infinity: 

lim sup \{p — p,u — u,9 — 9)(t,x)\ = 0. (1-24) 

To derive the large-time behavior of solutions to the compressible Navier-Stokes equations (1.1), 
it suffices to deduce certain uniformly-in-time a priori estimates on the perturbations toward the 
asymptotic state and the essential step is to obtain the positive lower and upper bounds on the 
density p{t, x) and the temperature 9(t, x) uniformly in time t and space x. In the case of small 
perturbation, one can use the smallness of the a priori iJ^-norm of the perturbation to get the 
uniform bounds of the density p and the temperature 9. Owing to such uniform bounds and the 
smallness of the boundary strength 5, one can derive certain uniform a priori energy-type estimates 
as shown in [15, 26]. In the case that the adiabatic exponent 7 is close to 1, by observing that 
9 = for ideal polytropic gases (1.I)-(I.2), one can deduce that \\9 — 1||lo°([o,t]xR) 

can be sufficiently small. Thus the desired energy-type a priori estimates can be performed as in 
[25, 26] based on the smallness of 5 and the a priori assumption 

\ < 9{t,x) < 2 for all {t,x) G [0,T] x M. 

However, these arguments are no longer valid for the case with large initial perturbation and 
general adiabatic exponent. We note that, even for the asymptotic stability of constant state to 
the Cauchy problem for the system (I.l), the uniform positive lower and upper bounds on 9{t,x) 
are given only very recently by Li and Liang [18], although the corresponding uniform bounds on 
p{t,x) were addressed in [12, 13] thirteen years ago. In their work [18], Li-Liang considered the 
fixed-domain problems to the compressible Navier-Stokes equations in the Lagrangian coordinate 
and obtained the uniform positive lower and upper bounds on the temperature 9{t,x) through 
a time-asymptotically nonlinear stability analysis. However, the outflow problem (I.l)-(1.5) will 
be transformed into a free boundary problem in the Lagrangian coordinate, which makes the 
treatment of boundary more difficult. To overcome this difficulty, we shall make use of a direct 
energy method to to the reformulated problem for the compressible Navier-Stokes equations (I.l) 
in the Eulerian coordinate and take account of the dissipative effect of the boundary terms. 

The main point for deriving our main results, the stability of the rarefaction wave and its 
superposition with a non-degenerate stationary solution to the initial boundary value problem 
(1.1)-(1.5), is to employ the smallness of the boundary strength 6 to control the possible growth 
of the perturbation suitably. Specifically, we first deduce the basic energy estimate with the aid of 
the decay properties of the smoothed rarefaction wave and the non-degenerate stationary solution 
provided that the boundary strength 6 multiplied with a certain function of mi (the a priori 
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lower bound of density p), m 2 (the a priori lower bound of temperature 9) and N (the a priori 
bound of the L°°{0, T ; iJ^(IR_|_))-norm of perturbation) is suitably small (see Lemma 3.1 for detailed 
statement). Next, to get uniform pointwise bounds of the density p(t, x), we transform the outflow 
problem (1.1)-(1.5) into a free boundary problem in the Lagrangian coordinate and modify Jiang’s 
argument for hxed domains in [12, 13]. Especially, we will use a cut-off function with parameter to 
localize the free boundary problem, and then we will deduce a local representation of the specihc 
volume V = 1/p to establish the uniform bounds of v. With such uniform bounds of the density 
p in hand, we can derive the iJ^-norm (in the spatial variable x) estimate of the perturbation 
uniformly in the time t in the Eulerian coordinate. And the maximum principle enables us to 
get the positive lower bound of the temperature 9{t, x) locally in time t. In view of the a priori 
assumption (3.8), we have to obtain the uniform positive lower bound of the temperature 9{t,x), 
which will be achieved by combining the locally-in-time lower bound of 6(t,x) and a well-designed 
continuation argument. 

The layout of this paper is as follows. After stating the notations, we summarize the existence 
of the stationary solution and some properties of the smoothed rarefaction wave in Section 2. The 
basic energy estimate, the uniform bounds of the density p, the uniform iJ^-norm estimate and 
the locally-in-time lower bound of the temperature 9 will be obtained in subsections 3.1, 3.2, 3.3 
and 3.4, respectively. The last part of this manuscript, subsection 3.5, is devoted to showing the 
proof of our main results by applying a well-designed continuation argument. 

Notations. Throughout this paper, L'^(IR+) (!<(?< 00 ) stands for the usual Lebesgue space on K_|_ 
equipped with the norm ll-llLq and iL^(K_|_) {k G N) the usual Sobolev space in the sense with 
norm || • ||fe. We introduce || • || = || • for simplicity. The space of continuous functions on the 

interval / with values in iJ^(]R+) is denoted by C{I] or simply by C{I; H*^) while the space 

of L^-functions on I with values in iL^(]R_|_) is denoted by L^(/; i?^(M+)) or simply by L^(/; H^). 
The Gaussian bracket [x] means the largest integer not greater than x, and x+ := max{a;, 0} is the 
positive part of x. 

2 Preliminaries 

It is well-known (see [2, 29] for example) that for each {t, x) G R+ x K, 

{p^,u^,9^){t,x) G R^[p-,u-,9-), X3{p^,u^,9^){t,x) = w^{t,x), 

where 

A 3 (p, u, 9) :=u+ \/Rr/9 

is the 3-characteristic speed of the system (1.9) and w^[t,x) is the continuous weak solution of 
the Riemann problem on Burgers equation 

J wf' + = 0 for (t, x) G M+ X R, 

\ w^{0,x) = w± for ± a; > 0 

with w± = X3 {p±,u±,9±). Moreover, w^{t,x) takes the form of 


' W- 

for x < W-t, 

x/t 

for w-t < X < w+t. 

^ w+ 

for X > w+t. 


The main idea in [8] is to approximate w^it, x) by the solution wit, x) of the Cauchy problem (1.11). 
The following lemma can be deduced by virtue of the method of characteristics (see [8, 23]). 
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Lemma 2.1. Let W- < w+. Then the Burgers equation (1.11) has a unique smooth solution 
w{t,x) satisfying 

(i) Wx{t,x) > 0, W- < w{t,x) < w+ for (t,x) G M+ x R; 

(ii) when x < W-t, w{t,x) — W- = Wx{t,x) = Wxx{t,x) = 0; 

(iii) for each p G [l,oo], there exists a eonstant Cp^q sueh that 

lka:(i)|lLP(R) < Cp^qmm{w,wh~^~''p}, 

(iv) lim 4 _>oo sup^gR \w{t, x) - w^{t, x)\ = 0 . 

Having obtained w(t,x), we can define the smoothed rarefaction wave {p,u,9) according to 
(1.12). Then one can check from a direct calculation that {p,u,9) solves the compressible Euler 
system (1.9) and 

(p, It, 0)(O, t) = (p_, M_, 0_), lim (p, It, 0)(t, x) = (p+, 0+) for each t > 0. (2.1) 

x—>-oo 

In view of (1.12) and Lemma 2.1, we have the following properties for the smoothed rarefaction 
wave (p, u, 9). 

Lemma 2.2. The smooth approximation {p^u,9) connecting {p-,U-,9-) and (p+,u+,0-|-) satis¬ 
fies 

(i) p_ < p < P+, 0_ < 0 < 0+, 9x(t,x) > 0, and 

P. = -^p9-^x, Ux = ^^0-^0.; (2.2) 

7 — 1 7 — 1 

(ii) if X < {u- + ^/Brf9Z){l + t), then {p,u,9){t,x) = {p-,U-,9-); 

(iii) for each p G [1, oo], there exists a constant Cp^q such that 

||(Pa:,ax,0x)(i)7p(R^) < C'p,,min{(5,57l + i)”^’^^}, (2.3) 

\\{Pxx,Uxx,9xx)it)\\p^p^^_^^ < Cp,, min{(5, (1+ (2.4) 

where S := |(m+ — U- , 0+-0 _)| is the boundary strength; 

(iv) lim 4 ^ooSup,,gR^ |(p,M,0)(t,x) - (p^,u^,0 ^)(l+t,x)| =0. 

Next we state the existence and the properties of the stationary solution (p, ft, 0) satisfying 
(1.16) and (1.17), which has been derived in [15]. To this end, we introduce the Mach number at 
infinity as 



Cm 


where Cm '■= \/Rl9m is the sound speed. 

Lemma 2.3 ([16]). Suppose that (u_,0_) satisfies 

(u_,0_) eTW™ := {(u,0) 


for a certain positive constant (5o. 
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(i) For the case > 1, there exists a unique smooth solution {p,u,9) to the problem (1.16)- 
(1.17) satisfying the decay estimate (1.18). 

(ii) For the case Mm = 1, there exists a certain region C Al™ such that if (u-,9-) € 
then there exists a unique smooth solution {p,u,9) to (1.16)-(1.17) satisfying 

f<Xn+l 

\d'f{p- PmW- Um,9 - 9m){x)\ < - - - + CSe~^^ for all n £ N. 

(1 -I- dx)'^+^ 

(ill) For the case Mm < 1, there exists a curve A4~ C Al™ such that if {u-,9-) £ Al"*, then 
there exists a unique smooth solution {p,u,9) to the problem (1.16)-(1.17) satisfying (1.18). 

Since the stationary solution (p, u, 9) and the smoothed rarefaction wave (p, u, 9) are well- 


defined, one can deduce that (p, u, 9) satisfies 

{ Pt + upa, -I- pux = A, 

p{ut + UUx) + Px = pUxx + f 2 , ( 2 - 5 ) 

Ovp{9t -\~ u9x') Pux — f^9xx A h^x A A 

for (t, x) £ 1R+ X K and the condition 

{p,u,9){t,0) = {p-,U-,9-), lim {p,u,9){t,x) = {p+,u+, 9+) for each t > 0. (2.6) 

x—¥oo 

Here P := P{p, 9) = Rp9 and 

fl = Ux{p - Pm) A Uxip - Pm) A pxiu - Um) A Pxiu - Um), (2.7) 

/2 = p[Ux{u - Um) + Uxiu - Um)] + UUx{p - Pm) 

+ (2.8) 
/a = Cyp[9x{u - Um) A 9xiu - Um)] A Cy{p- Pm)u9x 

A Ux{P - P) + {P - P)Ux - R9[p - Pm)Ux- (2.9) 


3 Stability analysis 

This section is devoted to proving our main results: Theorem 1 and Theorem 2. We will 
concentrate on the proof of Theorem 2, that is, the stability of the composition of a rarefaction 
wave and a non-degenerate stationary solution. The proof of Theorem 1 is similar to and simpler 
than that of Theorem 2. We therefore omit it here for brevity. 

First we introduce the perturbation {(f, if, d) toward the superposition wave (p, u, 9) as 

{(j),ilj,'d){t,x) := {p,u,9){t,x) - {p,u,9){t,x), 

where {p,u,9) is given by (1-21). Then we subtract (2.5)-(2.6) from (1.1)-(1.5) to have the initial 
boundary value problem: 

{ 0t A u(j)x A pifx = /i, 
p('0t A UTpx) + {P - P)x = ppxx A /2, 

Cypi'dt A wdx) A Pifx = K-dxx A pifl A h 


(3.1) 
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for {t, x) € K.+ X R+ with the initial and boundary conditions 

= (0,0). (3.2) 

Here the initial condition (<^ 0 , V'o, ^?o) := (po,uo,^o) ~ (/3,u, 0)|t=o satisfies 

lim (())o, V'o, = (0, 0,0), (3.3) 

X^OO 

and 

fi= - Ux(j) - PxiJ - A, 

/2 = PUxx - p.p~^(l)Uxx + P~^<t>Px - pi^Ux - P~^pf2, 

/a = - p~^4>{k0x;x + pul) + pul 

+ 2p'>pj:Ux: + 2pU^Ux - Rp0U^ - Cyp0j:1p - P~^pfs, 

where A (* = 1,2,3) are dehned by (2.7)-(2.9), respectively. 

We turn to deduce some desired a priori estimates for the perturbation (i^, ip, d) in the Sobolev 
space . Before doing so, for some non-negative constants N, s, t and rm {i = 1,2) with t > s, 
we introduce the set in which we seek the solution of the initial boundary value problem (3.1)-(3.2) 
as follows 

X{s,t;mi,m 2 ,N) := {{(p,ip,'d) e C{[s,t];H^) : (ipx.dx) G L^{s,t-,H^),(px, G L‘^{s,t-,Lp-), 

\\{(p,ip,d){t)\\i < N, {(p-\- p){t,x) > mi, {d-\-0){t,x) > m 2 V(t,a:) G [s,t] x M+}. 

The letter C or Q {i G N) will be employed to denote some positive constant which depends only 
on inf 2 ;gR_,_ {po(a^), ^ 0 ( 2 ;)} and ||(0o, V'o,'*^o)l|i- The exact value denoted by C or Ci may therefore 
vary from line to line. For notational simplicity, we introduce A ^ B '\1 A < CB holds uniformly 
for some constant C. The notation A B means that both A < B and B < A. Besides, we will 
use the notation (p, 0) = {cp p,!) 0). 

To make the presentation clearly, we divide this section into five parts. The first four parts 
concern the a priori estimates for the solution {(p,ip,d) G X(0,T;mi,m2,iV) to the problem (3.1)- 
(3.2), where T > 0 and it will be assumed that mi < 1 < N (i = 1, 2) so that 

\\{(p,ip,'0)(t)\\i < N, mi < p{t,x) < IV, 7712 < 0(t,x) < N for all {t,x) G [0,T] x R+. (3.7) 

In Subsection 3.5, the last part of this section, we will combine the energy estimates with a well- 
designed continuation argument to prove Theorem 2. 


(3.4) 

(3.5) 

(3.6) 


3.1 Basic energy estimate 


In this part, we will show the following basic energy estimate. 


Lemma 3.1. Suppose that the conditions listed in Theorem 2 hold, 
small Co > 0 such that if 

E{mi,m2,N){S -\-5) < cq 
with S(mi,m 2 ,iV) := then 


Then there exists a sufficiently 


(3.8) 


sup 

0<t<T 


tig 

0 0 ^ 


^ + -Ar da;dt < 1, 


f pSdxf p$ (t, 0)dt-I- f f 
Jk.+ Jo \pJ Jo J'r 

sup f ^dx+ r ^{t,o)dt+ r [ ^dxdt<N, 

0<t<TjR, P'^ Jo P Jo Jm^ P 


(3.9) 

(3.10) 
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where 


£ := ^ 




$(z) := z-lnz- 1 . 
Proof. Step 1. After a straightforward calculation, one can derive 


{p£)t 

from which we obtain 




pu£ + tjj^P - P) - mV’V'x - 


9=1 


^ / p£dx-u-p^(^\{t,Q)+ [ 

QWr, VP/ /Rh 




01?? 




dx = j2j 

9=1 "'R+ 


T^-^dx, 


(3.11) 


(3.12) 


(3.13) 


where each term TZq on the right-hand side of (3.12) will be defined below. Before defining and 
estimating all the terms on the right-hand side of (3.13), we set 

U . (p, 11 , 6 '), tj . (p, U, 0), Um ■ iPm: ; 0m)? ^ • (^? 1/?? 0)- 

First we consider 


7ei := 


0a0a , , 

K—-h 2pUxiPx 


which is trivially estimated by Sobolev’s inequality as 

\T^i\<mf^\Ux\\m'^x\<mf^d!\\^\\d>x\\^\Ux\\d!x\. 

In light of (2.3) and (3.7), we apply Holder’s and Young’s inequalities to deduce 

[ |7^l|dx<m2-2||lI/||^|lC/a||||lE'a||^ 

/r+ 

< -I-1)“^ ||iE'a;||^ 

< {l + t)-^ + m~^Nisi\\d'x\\^. 

Next we consider the term 

0 

77-2 .— ~r{i^6xx d~ pa^') P pipUxx- 
a 

This term can be controlled by Sobolev’s inequality as 

\T^ 2 \<mf^\d>\ [\Uxx\ + \Ux\'^] [\Uxx\ + \Ux\^] 

According to (2.3), (2.4) and (3.7), we infer 

[ |7^2|dx < m^-^llibll^lliball^ [WUxxWl^ + \\Uxf] 

Jr+ 

< to^^N^||iI'2;|1^0^(1 -I- t)~i 


(3.14) 


(3.15) 


< 


(1 + t) 6+7712 N II 47 


Let us now consider the term 
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It is not difficult to derive from (2.7)-(2.9) that 

lifukM < m\u - u^\ + \k\\u - Uml 

In view of Lemma 2.2, we deduce that U{t, 0) = Um and hence we have that for g > 1, 

\k\\U - Ural + \U.\\U - Um\ + \U^\k 

< 


Li 


\U.\\U -Ura\ + \U,\ / \U^\{-,y)dy + \U,\\U, 

do 


Li 


LI 


<||C/x||l» \U -Urr.\+X\k\ + \U. 

which combined with (1.18) implies 

\UA\U - Ural + WAIU - Ural + W^lUxl < 

Li 

It follows from (3.7), (3.16) and (3.17) with g = 1 that 

[ |7^3|da:<iVm^^||^'|Uoo5||i/,||io.. 


(3.16) 


(3.17) 


We use (2.3) and Young’s inequality to have 


|7^3|cl:r < Nm^^Sil + t)-id /||2 ||vl/,|| 2 < (l + t)-3 + . 


(3.18) 


We then estimate the term 


77.4 := tjj9a; 


R(j)-\- Rp^ ( y ) + ( 3 


0 


PUx 


+ R9^ 

Cv \pj \9 


pi’Ux 


1 - 


+ i7 


1 _ _P 

p9\ L' 


P 

P L 


k9xx + pul 
n9xx + pul + fz 




9x'9x 

K— - 1 - 2pUxll^x 


(j) 


To this end, we first obtain from the identity 

r-l /•! 


$(z)= [ f 9i^"il + 9i92{z-l))d92d9i{z-lf 
Jo Jo 


that 


^-2/ 


^-1 


This last inequality implies 

P 
P 



(3.19) 


(3.20) 
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In light of (3.16) and (3.20), we discover 


[ |7^4|dx < wee / \di{U-Urn)H^\^dx. (3.21) 

Jr+ ^^g JR+ 

To estimate the terms on the right-hand side of (3.21), we utilize an idea in Nikkuni-Kawashima 
[24], that is, the following Poincare type inequality 

\ip(t,x)\ < |<p(t, 0)1 -I- y/x\\ipx{t)\\ for X e M+. 

Applying this inequality to 41, we deduce from (1.18) and (3.19) that 


\diiu - Um)\md>\^dx < d</>(t,0)^ + S\\^4^ < N^m^^Sp<h ^ (t,0) + S\\^,r (3.22) 


for fc = 0,1 and £ G N. Plug (3.22) into (3.21) to deduce 


|7^4|da; < 


We next estimate the term 


Tib ■= i’Ox 


-—pd + R(l)+ Rp^ ( - ) + CvP^ ( - 


PUx 




$ ( - ) + ^ + RSd> 


To this end, we introduce 


In light of (2.2), one can Hnd 


a := In - and 6 := In ( - . 


Tib = -p0xF[ip,a,b) 


(3.23) 


(3.24) 

(3.25) 


with 

a, b) := Rij^a + Cvi^b + \JR'jR0~^0{e°' — a — 1) 

+ ^ y^cj-i9{e^ -b-1). 

7- 1 

One can easily deduce that {F, d^F, daF, dbF){0, 0,0) = (0,0,0, 0) and that the Hessian matrix of 
F at point (0, 0,0) is 


^2 9 2 i? Cy 

Hf= R y/WjR9-^9 0 

\ Cy 0 ^JR'ycy9~^9/ 

Let Dk (A: = 1, 2, 3) be the k x k leading principal minor oi Hp. Then we have 


D.=2kS»-4, D, = cV^R‘i-i W^"‘'-\ 

7 — 1 7 — 1 7 — 1 
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Apply Sylvester’s criterion (see [5, Theorem 7.2.5]) to deduce that there exists a positive constant 
El such that if (5 < ei, then Hp is positive definite. Then it follows from the Taylor formula (see 
[1, Theorem VII.5.8]) that 

2 1 

F{'tp,a,b) = ^ ada-\-bdb)’^F{0,0,0) 


9=0 


1 

+ - {1 - t)'^{^|^^^ada + b^b)^F{t^j;,ta,tb)dt 

^ Jo 

< l-{'tlj,a,b)F[p{'ip,a,b)^ + f (i?a^e*“ + c„6^e**’) dt. 

2 2 Jq 


Due to 9x > 0, we derive that if <5 < ei, then 

7^5 < PO: 

We note that the identity 


Inz = 


a2(ela| _ i) + 52(el^l - 1) 
(z - I)d6li 


/o 


1 F 0i{z — 1) 


implies 


|lnzp< (z-i + l)2(z-I)2. 

Then we apply (3.27) to a and b and use the estimate (3.26) to find 

which combined with (2.3) yields 

/ |7^5|dx < 

< (1 + t)-i II ^111II ^,11 i 

< (l + t)“® + ■ 

Plug (3.14), (3.15), (3.18), (3.23) and (3.28) into (3.13) to obtain 


^ / p£dx-u_p^(^\{t,0)+ ( 
dIjR, \pj Jr 


^ + d. 
0 ^ 02 


< (1 + t)-s + 


~5^ +5 


p4> 




0) + ll'I'x|P 


Hence we can hnd a sufficiently small constant £2 > 0 such that if 




+S 


< £2, 


then 


_d 

dt 


pSdx — pu-^ ( - I (f, 0) + 


< 


(1 + t)"s 


-'R^. 


ib'^ 02 
zA + da; 
0 02 


U4\ 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


(3.30) 
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from which we obtain 


sup / p£dx + 


0<t<T. 


r -0 


{t, 0)dt - 


6» 6»2 


1 I —12 —12 7vrl6 

1 + mi 7712 


(52 + d 


Jo 


(3.31) 


dt. 


Step 2. We now make some estimates for the last term in (3.31). We first differentiate (3.1)i with 
respect to x and then multiply the resulting equation by (j)x /to find 


2p3 


'^4^'x \ ~ Jx , ^ (j^x'^px _ „ 4>x (px'p’x 

r, 3 I ~ 3 — 2~ 

V/x P P P P 


Multiply ( 3 . 1)2 by (px !to have 


( Px1p\ 

PtJ’ Px'tp'^ 

\ p Jt 

. p p . 


P^xP^xx 


PxUxpIp Pxx ,2 Pxfp'px Pxlpfl PxfppJx 


- ^ _ 

p p p 


Ux'i’xP’ 'P’x t ,2 I Uxpx'p , JJ £,1 <Px r (px 

- + -/l +'P’x + - {P-P)x—+ /2 —. 

p p p p p 


In light of (3.32) and (3.33), we have 


(3.32) 


(3.33) 


' ppl , pxlp' 


' pupl 

Pti’ Px'P'^ 

1 

f 

Iro 

1 

_ 1 

+ 

t 

2p^ 

1 

1 

1 


Repl 


= E«.- 


q=l 


(3.34) 


which combined with Cauchy’s inequality implies 
rt 


f r4(s,0)ds+ r / M<1+/ pP;^dx + J2f[ I2?l. (3-35) 

dR+ P Jo P Jo Jr+ P Jr+ ^iJo Jr+ 


where each term Qg on the right-hand side of (3.34) will be defined below. First, let us define 


Qi := rx - 


2 R4^x'^x 


and by applying Cauchy’s inequality, we have 


/7 iq.i<77 21? 

Jo J^^ Jo J^ 


t r V .q2 


(3.36) 


Then we consider the term 


Q2 := - {uxxp + Pxxip) - ^^{puxxp - Rpxd(p - UxPp'tp) 

(p4>,~ . , - - N Pxxi’’^ ^Pxlplpx , Uxplpx Rpxdpx 

- :r[PxUx + PxUx) - 2 -1- 5 - , 

P P P P P 
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which can be controlled as 

I22I < rn-^{S + 
In light of (3.22), we have 


For the term 


\Q 2 \^N^m^ S pd) - (s,0)ds + TO^ (i5 + (5) / ||vl/ 3 ,(s)|| ds. 


Qs ■■= - ^^^{Uxx(t> + Pxxtp) - ^^{-RpxH - UxP(t>ll^) 




p^p 


- 2 


Pxi’il’x , Ux4>i^x Rpxd4>x 


p 


p 


p^ 


it follows that 

\Q3\<m^''\{Ux,Uxx)\ l^ll^.l \{Ux,U..)\ l^x| • 

In view of (2.3), (2.4) and (3.7), we apply Holder’s and Young’s inequalities to deduce 


|Q3l<m-3 / ||vl/|l^|l(i7,,f/,,)||||vl/,||ids 

pt 

/ (1 + = ||4'a;(s)|| 2 ds 

Jo 

< 1f ||4'2 ;(s)||^ ds. 

Jo 


For the term 


we have 


24 :=- 


PxUx(j)'lp Pxxi’'^ 


I 24 I < mr'|4/p [\Uxx\ + \Ux\^] < mr'II 4 /IIII 4 /,II [|C7,,| + |C7,|2] , 
which combined with (2.3)-(2.4) and (3.7) yields 


r / |24| < mr" f Ill'llll'E'.ll [WUxxWl^ + \\Uxf] ds 
Jo Jr^ Jo 


<mi^N / ||4'a,(s)||5«(1 + s) '‘ds 


<l + m-{^N‘^5^ I ||4'^(s)fds. 


Finally we consider the term 


^ p4>x ; 4‘x , _ ^ , Px'i’fi ''Pxfi 

25 :=- ^Jix - ^(-puxxp + PJ2) - 2 -^-■ 

P^ P P P P 


Holder’s inequality gives 


rt r 


\Qs\<ni^^ / II'I^.IIII(«..,/i,/2,/i.)II + II«'IIl<»II/iIIl‘ 


10 


ds. 


(3.37) 

(3.38) 


(3.39) 


(3.40) 


(3.41) 











16 


L. Wan, T. Wang and H. Zhao 


(3.42) 


\U,\\U^\ + |C7,,||t7 - U^\ + |C7,,||C/ - U^\ < (5||(C7„ t7,,)||ioo. 

LI 


We deduce from (2.7) that 

\fl.\ < \U.\\U.\ + \U,,,\\U - Um\ + \U..\\U - Urr 
Similar to the derivation of (3.17), we have that for q > 1, 

(3.43) 

Combining the estimates (3.16), (3.17), (3.42) and (3.43) and utilizing Lemma 2.2, we deduce 

11 id^xx 1 f 1 ^ f2: f Ix^ ^ II '^xx ('^) 11 11 i^x : U^x^ ('^)||j^oo ^ ^’^(1 “t“ s) 2®. 

Plug this last estimate into (3.41) and use (3.16)-(3.17) again to have 

IQsl < 

< 1 + 


(mi 


|vI/,pr(l + s)-^+<5(l + s)-l||vI,||2||vI;,^||5 

rt 


ds 


+m7^^N^~6^ 


(3.44) 


|4'x(s)|pds. 


Plugging (3.31), (3.36), (3.38), (3.39), (3.40) and (3.44) into (3.35), we take e sufficiently small to 
find 






< 1 + 


oP^ Jo J«.+ P 


s,0)ds+ f f 
Jo JR^ 


rx + 


(5= + (5 


- 12 .^- 12^16 


+ rrii m2 




Il<(’x(s)f ds 




which combined with (3.31) gives 

/■* ^2 
|dx+ / 

K+ P 


dx+ [ 4(s.0)ds+ [ / 

Jo P Jo Jr 


0 Jr+ P 


< 1 + 


0 JR+ 


€ 


9 




5^+5 


0 Jr+ P 


We take cq > 0 small enough and use (3.8) to have 


Jr+ P Jo P Jo Jrx P Jo Jr 


0 Jr+ P 


10 ^R+ 


q2i 


rx + f 


(3.45) 


The estimate (3.9) follows by plugging (3.45) into (3.31) and using the condition (3.8) for a suf¬ 
ficiently small eo > 0. Combine (3.9) and (3.45) to deduce (3.10). The proof of the lemma is 
completed. □ 


3.2 Uniform bounds on density 

Having obtained the energy estimate (3.9), we can proceed to deduce the positive lower and 
upper bounds of the density p{t,x) uniformly in time t and space x in this subsection. For this 
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purpose, we transform the outflow problem into the corresponding problem in the Lagrangian 
coordinate by introducing the Lagrangian variable 


y = -u- 


[ P{s,0)d 

Jo 


s+ / p{t,z)dz. 

'o 


(3.46) 


By the coordinate change {t,x) i->- (t, y), the domain [0,T] x R+ is mapped into 

rtr ■= {{t,y) ■ 0 di t < T,y > Y (t)} with Y(t) := —u- f p(s, 0)ds, 

Jo 

and the outflow problem (1.1)-(1.5) is transformed into the following initial boundary value problem 


Vt-Uy = 0, 

U.dPy=(^) , 

\ V Jy 

. + + for!,>y((), (3-47) 

(a, ^)|y—v(i) 

{v,u,e)\t^a = {vq.uo.Oq). 

Here v = 1!p stands for the specihc volume of the gas and vq = 1/po- The basic energy estimate 
(3.9) in Eulerian coordinate can be rewritten as a corresponding estimate in Lagrangian coordinate 
as a direct consequence of the transformation (3.46). 

Corollary 3.2. Suppose that the conditions listed in Lemma 3.1 hold. Then 


sup / £dy + 
o<t<T Jvit) 


pT poo 


Jo JYit) 

1 

bO 


dydt < 1. 


(3.48) 


Note that the function Y {t) describing the boundary in the Lagrangian coordinate is part of the 
unknown, that is, the problem (3.47) is a free boundary problem. To obtain the uniform bounds 
of the specific volume v for the free boundary problem (3.47), we introduce the time-dependent 
domain Hi{t) with i € Z and t G [0, T] as 


Lli{t) 


[Y{t),[Y{t)]+2] dz = [Y{t)] + l, 
else. 


(3.49) 


Based on the basic energy estimate (3.48), we have the following lemma. 

Lemma 3.3. Suppose that the conditions listed in Lemma 3.1 hold. Then there exists a positive 
constant Co, depending solely on infa;gR^{/3o(a;),0o(®)} cind ||((^o,'0o,^?o)||i, such that for all pair 
(s, t) with 0 < s < t < T and integer i > [Y (t)] -|- 1, 

< f v{s,y)dy <Co, < f 6»(s, y)dy < Co, (3.50) 

and there are points ai{s,t),bi{s,t) G Tlift) satisfying 

Co"^ ^ v{s,ai{s,t)) < Co, Cq"^ < 8{s,b^{s,t)) < Cq. 


(3.51) 
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Proof. Let 0 < s < t < T and i > + 1- According to the definition of Y{t) and the sign of 

u_, we have L"(s) < Y{t) and C [F(s),oo). In view of (3.48), we get 

[ ^ (t) (s,y)dy + / $ (s,?/)dy < 1. 

Apply Jensen’s inequality to the convex function $ to obtain 


$ 




-{s,y)dy 

V 


+ 4* 



iniit) 0 


y)dy 


< C. 


Let a and j3 be the two positive roots of the equation $( 2 ) = C. Then we have 


a < 


1 

WW\ 



V 

-{s,y)dy < P, 

V 




These estimates imply (3.50). Finally we employ the mean value theorem to (3.50) to find 
ai{s,t),bi{s,t) G satisfying (3.51). The proof of the lemma is completed. □ 


We deduce a local representation of the solution v for the free boundary problem (3.47) in the 
next lemma by modifying Jiang’s argument for fixed domains in [12, 13]. To this end, we introduce 
the cutoff function G IF^’°°(R) with parameter 2 G R by 


1, y < [ 2 ] +4, 

Pziy) = {[z] + 5-y, [ 2 ] + 4 < y < [ 2 ] + 5, 

0, y> [ 2 ] + 5. 


Lemma 3.4. Let {t,z) G fir- Then 


v{t,y) = B„{t,y)A„{t) 


R B,{t,y)A,{t) 


y Jo B„{s,y)A„{s) 
for all t G [0,t] and y G Iz{t) := (Y(t),oo) fl ([2] — 1, [2] +4), where 


9{s,y)ds 


Bz{t,y) :='!;o(y)exp|iy (mo( 0 - ^(t,0) 7’z(0dc| , 

A.(t):=expjir r\^-p)dfds]. 

[yJo J[z]+4 ^ V J j 


Proof. We multiply (3.47)2 by ipz to get 


-y,'AyJ^-P). 

l\ V \ V / 


(3.52) 


(3.53) 


(3.54) 

(3.55) 


(3.56) 


Let {t,y) G [0,r] x Iz{t). Since y > l^(s) for each s G [0,t], we have [0,t] x [y,oo) C IIt- In light 
of the identity ipziy) = 1 and (3.47)i, we integrate (3.56) over [0,t] x [y, 00 ) to get 
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This implies that for each t G [0,t], 

1 ^ fR \ 1 

i’{s,y) B^{t,y)A^{t)' 


(3.57) 


Multiplying (3.57) by R9{t,y)/fi and integrating the resulting identity over [0,t], we have 


exp 


h Jo v{s,y) J p Jo B^{s,y)AAs) 


We then plug this identity into (3.57) to obtain (3.53) and complete the proof of the lemma. 


□ 


The following lemma is devoted to showing the bounds of the specific volume u ( t , z) uniformly 
in the time r and the Lagrangian variable z. 


Lemma 3.5. Suppose that the conditions listed in Lemma 3.1 hold. Then 

< v(t,z) < Cl for all {t,z) G TIt. 

Proof. Let (r, z) G TIt be arbitrary but fixed. The proof is divided into three steps. 
Step 1. It follows from Cauchy’s inequality and (3.48) that 

B^{t,y) - 1 for all {t,y) G [0,t] x 7^(t). 

Let 0 < s < t < T. For each 0 < t' < t, there exists y{t') G [[ 2 ] + 4, [z] + 5] such that 

2/(^0) = ,,, inf 

(N+4,N+5) 


(3.58) 


(3.59) 


Apply Cauchy’s inequality to have 

p ^ R^ ^ I R<^ 

V ~ v9 2v ~ v9 v9 2v 

In view of (1.18), (2.3), (3.48), (3.8) and (3.50), we apply Jensen’s inequality for the convex function 
l/x to deduce 


t Azl+S 


fi 


s 7[z]+4 L "V 


pUy 


- P 


H r /-R+o 

< CCNmf^{SS){t - s) - — 9{t',y{t')) v~^{t',y)dydt' 

^ Js J[z]+A 

~ - B f* \ 7 W+® 

< CCNmf^{5d)(t — s) — — / 9{t',y{t')) / vdy 

2 Js [JW+4 

<C + Ceoit -s)- C-^ f 9{t', yit'))dt', 

J S 


(3.60) 


Since [z] + 4 > [F(t')] + 2, we derive from (3.49) that n[ 2 ]+ 4 (t') = [[ 2 ] + 4, [z] + 5]. We then apply 
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Holder’s and Cauchy’s inequalities to obtain from (1.18), Lemma 3.3 and (3.48) that 
i-t rvd') a 

' / (J.I 


If 

US b\, 


Is ^ 


< 








d^dt' 


< 7712 ^(^ + ^)(t — S) + 




' t^[2;]+4(t') 


v{t' 


dt' 


(3.61) 


/■* r°° 

<C{t-s) + C / 


Js JY{t') 

< C{t -s) + C. 

Applying Jensen’s inequality to the convex function e^, we have from (3.51) and (3.61) that 
f 0(,t\yit'))dt'= [ exp(ln6»(t',y(t')))dt' 


> (f — s) exp — J \n9{t',y{t'))dt' 

rvit') 


> [t — s) exp 


t — s 


/h[^j+4(i'.i') ^ 


^ (^', + In0(1', 6[,,]+4(t', f')) 


dt' 


> {t — s) exp — In Co — 


t-s f C 


t — s 


t pyit') a 

Is •t&[^]+4(t',t') ^ 


- f 0(t',y{t'))dt'< 

J S 


(3.62) 


j•ld+5 

_p 

Ad+4 

- V 


This implies 

0 if0<t — s<l, 

— C~^{t — s) if t — s > 1. 

Plugging (3.62) into (3.60) and taking eo > 0 small enough, we have for each s G [0, t] that 

< C-C-i(t-s). 


According to the definition (3.55), we then obtain 

0 < Az{t) < Ce”*^*^, < Ce“d-'’)/c fQj. all 0 < s < t < r. (3.63) 

Step 2. Plugging (3.59) and (3.63) into (3.53), we infer that for all {t,y) G [0,r] x /z(t), 

f 444^'(s,y)ds < 77(1,7/) < 1 + / 9{s,y)e~"~^ds. (3.64) 

Jo A^s) Jo 
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In light of the fundamental theorem of calculus, we deduce from (1.18) and (3.50) that for y G Iz{t) 
and 0 < s < t < r, 


6»(s,2/)2 - 6»(s,6[^]+2(s,t))2 


< 


e. 


(s,e)de 


< 7712 ^ (^ + '^) + 


< 7772 ^ (<5 + (5) + sup 77 2 (s, •) 

Iz{t) 




9 2|^j^|(s,^)d^ 


f 0?/ 

2 

f 



Jl.ir) 


(3.65) 


’Iz{r) 




where we have used h]^z]+ 2 {^i^) ^ ^ Combine (3.65) with (3.51) and (3.8) to give 

f f 'd^ 

1 - C sup 7;(s, •) / -^(s,0d'f ^ 0(5,7/) < 1 + sup 7;(s,-) / -^(s,C)dC. (3.66) 

We plug (3.66) into (3.64) to obtain 

j-t r ^2 

v{t,y)<l-\- sup7;(s, •)/ -^(s,0d^ds. 

Jo I At) J I,(r) v9 

Taking the supremum over /z(t) with respect to y, we have 

sup 7;(t, •) < 1 + / sup7;(s,-) / 

h(T) Jo I At) Jn 


'Qiir) 


^is,0d^ds. 


Applying Gronwall’s inequality to (3.67), we can deduce from (3.48) that 

sup v{t, •) < Cl for all t G [0, r], 

Iz{t) 


(3.67) 


(3.68) 


where Ci > 0 is some constant independent of t, r and z. Noting that z G Iz{t), we deduce from 
(3.68) that v{t,z) < Ci. Since (t,z) G fir is arbitrary, we conclude 


7(t, z) < Cl for all (r, z) G fir- 


(3.69) 


Step 3. On the other hand, in view of (3.50), (3.59) and (3.63), we integrate (3.53) on Izij) with 
respect to y to find 

ft 

1 < ' 


^ v{t,y)dy<e + f 444ds. 
lAr) Jo ^z{s) 


Consequently, we have 


Mt) 

Az(s) 


ds > 1 - Ce"*/^. 


(3.70) 
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Inserting (3.66), (3.69) and (3.70) into (3.64), we have 




Az{s) 

> 1 - - C 

> 1 - - C 


Azis) 

Jt/2 j ^z('S) 








f.t/2 


lO Jh{T) 

i-t r 


V 

V0^ 


dCds - C 


f I 

Jt/2 JIz{t) 


v9^ 


(3.71) 


^2 

vB^ 


> 1 - Ce"*/*^ - Ce-^ - C 

Jt/2Jh{r) 

> 1 for all (t,y) € [To,r] x 4(r), 

where Tq is a positive constant independent of t. In particular, the estimate (3.71) implies 

> 1 for all T >Tq, z > Y{t). (3.72) 

As in [16, 17], we can derive a positive lower bound for v, that is, 

v{T,z)>e~^* for (r, z) G Ht- (3.73) 

Finally, we combine (3.73), (3.69) and (3.72) to get (3.58). This completes the proof. □ 

As a corollary of Lemma 3.5, we obtain the bounds for the density p{t,x) uniformly in time t 
and space x. 

Corollary 3.6. Suppose that the eonditions listed in Lemma 3.1 hold. Then 

< p(t,x) < Cl for all {t,x) G [0,T] x IR+, (3.74) 

where the positive constant Ci depends solely on infa;gR^{po(a^),^*o(a;)} and ||((/)o, V’ 0 ;^?o)||i- 


3.3 Uniform estimates for the perturbation 

In this subsection, we will estimate the ijAnorm of the perturbation {(j),'d,'ip)(t,x) uniformly 
in time t. First we can get the following uniform L^-norm estimate. 

Lemma 3.7. Suppose that the eonditions listed in Lemma 3.1 hold. Then 

sup ||(((),'d,V^)(t)f + / [ [{l + e + tp‘^)tl}l + 'dl]dxdt<l. (3.75) 

0 <t<T Jo Jr+ 

Proof. We divide the proof into five steps. 

Step 1. First, for each t > 0 and a > 0, we denote 

:= {x G IR+ : 'd{t,x) > a}. 

Then it follows from (3.9) and (3.74) that 


sup 

0<t<T 



d^dx 





|(^(s,0)pds 



< C{a) sup 

0<t<T 



p£dx + / 
Jo 



(s, 0)ds 


+ C{a)f j 

Jo Js.+ 



^ X 
02 


< C{a). 


(3.76) 
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Step 2. We now estimate the integral f^, To this end, we multiply ( 3 . 1)3 by (i9 — 2)_|_ := 

max{il — 2 , 0 } and integrate the resulting identity over (0,t) x R+ to obtain 


f pit) - 2)\dxK f f dl-i-Rf f pOtjjxi'd - 2)^ 

Jrm Jo Jii2{s) Jo Jrm 

rt 


[ poido-2)ldx+ [ [ f,(d-2)++pf [ ^lid-2)+. 
JRa Jo Jr+ Jo JRa 


(3.77) 


0 ^R+ 


To estimate the last term of (3.77), we multiply ( 3 . 1)2 by 2'!/)(i? — 2)+ and integrate the resulting 
identity over (0,1) x IR+ to find 


[ 'ip‘^p{t}-2)+dx + 2p f f 'iPl{'d-2)+- f f pip^idt+ud^) 

JR+ Jo Jr+ Jo -'n^(s) 

= [ ijlpo{-do-2)+dx + 2R [ f peij4'&-2)++2R [ [ pO^d, 
Jr+ Jo v/r+ Jo Jci'^^s) 

+ 2R f [ i/jipeud - 2)+- 2p 
Jo JRa 


/O JR+ 

Combining (3.78) and (3.77), we have from ( 3 . 1)3 that 


i>tpxiJx + 2 / f2i>{iJ-2)+. 

0 ■In^(s) Jo Jr+ 


I 2 


p{d-2)l+i;^p{d-2)+ 


dx 


0 


+ M 


^"( 1 ^- 2 )- 


0 JR+ 


Cy 

12 


po(i?o - 2)++ V'qPo('!^o - 2)+ dx + ^Jp, 

p=i 


(3.78) 


(3.79) 


where each term J7p in the decomposition will be defined below. We now define and estimate all 
the terms in the decomposition. We first consider 


Ji ■■=R 


pBipxiiJ — 2)^ and J 2 '.= 2R 


/O ^R+ 

In light of (3.74) and (3.9), we have 


pOil^dx 


0 Jci'^is) 





pEdx < 1. 


(3.80) 


From Cauchy’s inequality and (3.74), we obtain 

\Ji\<ef [ ^l{d-2U+C{e) f [ e\d-2U 

Jo Jr+ Jo Jm.+ 

<ef j ^l[d-2)++C{e) f j 6{d-l)l 

Jo Jr+ Jo Jr+ 

<e/ [ 'Pli^-‘J)++C{e) [ sup(^»-l)+, 

Jo Jr+ Jo r+ 


(3.81) 
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and 


< e 


Here we have used the fact that 


\J 2 \<ef[ dl + C{e) f [ 

Jo Jo.2<,s) Jo Jfi'2{s) 

<ef j dl + C{e) f 

Jo J0,2(3) Jo J0,2(3) 

i f [ ^l + C{e) rsup(^-l)^. 
Jo Jn'(s) Jo R+ 


e < K{d - 1) 

with K = 2 + supjg^ 0 provided that d >2. Let us define 

rt 


Jo--= [ f [M^-2)+ + cZ^hi^^+2Mid-2)+]. 
Jo JnUs) 


/O Jfl'^is) 

According to (3.5) and (3.6), we use (3.16) and (3.74) to deduce 

l(/ 2 ,/ 3 )| <G+|(vi/,vi/,)||(i7,,t/,,c/,,)|. 

with 


G := \U^ + |t7,|" + \U,\\U - Um\ + \U.\\U - Um\ + \U,\\U.\. 


Hence, we have 


J3< 


G+|(4/,vI/,)||(c7,,t/,,;7,,)| [m-\d-2)l + ^^] . 


0 J 02 ( 3 ) 

It follows from Lemma 2.2 and (3.17) that 

||G'(s)||^i ^ (1 + s) 3 , 

from which we get 


ff G[m^\d-2)1 + ^^] 

Jo JOL(s) 


<S3 (1 + s)’ 




+ 7712^(53 / SUp(l? — 1) 
Jo R+ 


+ 7712 / SUp(7?— 1)^. 

Jo R+ 


Next we have from (3.7) and Lemma 2.2 that 

pt p 

|(4',^',)||(H„f7„C/,,)|(75-2)2 

|(C/.,H„t/,,)| \N^{d-2)l + \'i>, 


0 J 02 ( 3 ) 

pt 

-1 / 

mo 


)+ 


0 JQ.'2{s) 

<m^^N^{5 + S) [ [sup(79-l)^ + ||«', 
Jo R+ 


(3.82) 


(3.83) 

(3.84) 


(3.85) 


(3.86) 


(3.87) 


(3.88) 
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Since 


we have 


< (1+i^-i) 


Jo JO, 2 {s) 

<[ [ fiV|^|"|c7,|+iV2|^|'|(C7,,[/,,)|l +(^ + 5) / 

Jo Jqus) l -‘Jo 


< 1 + 


SbN 3 + (5iV^ 


/ + 5mr^iV4 / - (s,0)ds. 

/o Jo \Pj 


To derive the last inequality in (3.88), we have used (3.22), Young’s inequality and 

< Il'l'lll~rx|k4||vl/||l < 5^Ni(l + s)-™||'I/,||i. 
Plugging (3.87), (3.88) and (3.90) into (3.85), we deduce from (3.8)-(3.10) that 


|J3|<1+ [ sup(t5-1)2. 
Jo R+ 


Let us now consider the term 


Ji:=2R / ^PipeUd-2)+, 
Jo Js.+ 


which combined with (3.7) and (3.8) yields 


\J4\< I sup(^»-2)+||V'||||(p^, 0^)11 < / sup(^»-l)^. 


/O R+ 


/o M+ 


For the term 


J5 -= f f [hi’lcy - Pi^xc^ , 

J Q J fiofs) 


'0 ^ 02 ( 5 ) 

we apply Cauchy’s inequality and (3.80) to deduce 

pt p pt 


\Js\<ef j dl + Cie)f j V’V^+Z [ 

Jo Jfi'2{s) Jo Jn2{s) Jo dn^(s) 

<eff dl + Cie)f[ V’V^+r/ ^^(^-1) 

Jo da^(s) do dn^(s) Jo da^(s) 

< e / V € + C{e) f [ + f snp(^ - 1) 

Jo JviL(s) Jo JQt(s) Jo R+ 


)l- 


(3.89) 


(3.90) 


(3.91) 


(3.92) 


(3.93) 


Jo ■= [ [ c^^nii'^'dxx- 
J{) JQ.2 {s) 


We finally consider 
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In order to estimate we apply Lebesgue’s dominated convergence theorem to find 
Jg = — \im [ [ <p„{'d)'ip^da;x 

Cy V^0+ Jq 

= — lim / f [-2(p,,{'d)-ipil;x'&x - 
Cv i'^0+ Jo 

— lim / [ 2(p,,{-d)ilj'ipxdx 
Cv 1'^0+ Jo JR^ 


< - 


< e 


/ / € + C{e) 

JO Jr+ 


where the approximate scheme ^Pv{'&) is defined by 

('■ 

0 , 


d-2 >v, 

Q <d-2 <v, 

d - 2 < 0. 


Plugging (3.81)-(3.82), (3.91)-(3.94) into (3.79), we get from (3.74) that 

[ {d-2)ldx+ f f [dl + JjU^-2)+] 
jRj. Jo Jn'(s) 


< 1 + e 


f [ dl + C{e) f [ + C{e) f sup(^ - 1)^. 

do dR^ do dR+ do R+ 


0 K+ 


Step 3. We obtain from (3.9) that 


[ [ + j f [dl + i:l{'d-2)+]+ f [ 

Jo dR.,. do dn^(s) do dR^\n^(s) _ 

<f[ [dl + J;l{d-2)+]+l. 

Jo JCl'^(s) 

Combining (3.96) and (3.95), and choosing e sufficiently small, we have 


9^ e 


(3.94) 


(3.95) 


(3.96) 


/ {d-2)ldx+ [ [ [dl + J;l9]<l+ [ snp{d-l)l+ [ [ (3.97) 

Js.^. do dRo. do R+ do dRi 


Step 4. To estimate the last term of (3.97), we multiply (8.1)2 by and then integrate the 
resnlting identity over (0,t) x R+ to have 


ptp'^dx + 3/x 




1 


Po^oda; 


/o 


= 3R 


J) iJxQ4> + 3i? 

.+ do 


From (3.76) and (3.80), we have 


fj 

Jo Jr. 


Ip'^lljxOcj)+ 




Jj'^'iJxP'd 


(3.98) 


f2r 


0 dR+\a^(s) 


< 


[ll<^ll + IMU=(R+\n'(.))] < 

Jo Jo 


(3.99) 
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We then apply Cauchy’s inequality to derive 


0 


<ef [ + C'(e) f f 

Jo Jb+ Jo Jci'^^s) 

rt 


< € 


f [ + C'(e) / sup {d - 1)^ . 

Jo Jrj- Jo R+ 


(3.100) 


10 -'R+ -'O R+ 

In view of (3.83), we utilize (3.86), (3.89), (3.22) and the Young’s inequalities to have 


< N 

< 1 + 


G+|(«',^,)||(t7„t7„C/,,)| V’" 


+ 

-2 , ^ 34 


S9N- +SN^ 


^xlr + Sm^^N^ I ( ^ ) (s,0)ds. 


(3.101) 


Plugging (3.99) -(3.101) into (3.98), and taking e sufficiently small, we derive from (3.8)-(3.10) 
that ^ ^ ^ 

tj}'^dx-\- [ f [ -^1+ f sup(i?-l)+. (3.102) 

Jo JRj. Jo Jr^. Jo R+ 


It follows from (3.9) that 


f f f [ 0^1 + C{e) f f ^ < e r / + G(e). (3.103) 

JO JR+ Jo Js.4- Jo JRx " Jo JRa 


Jo JR+ 


JO JR+ 


0 JR+ 


Combination of (3.103) and (3.102) yields 


[ V’^dx-k / f il+'ip^)'ipl<C{e) + e f f 0 jIjI+ [ sup(i?-l)+. (3.104) 

jRj. JO JRj- JO JRi JO R+ 


JO JR+ 

We plug (3.104) into (3.97) and choose e suitable small to find 


{d - 2)1 + ^^dx+ / [dl + jl;l{l + 0 + ^2)] < 1 + / sup (zl - 1) 

Jo R+ 


0 JR+ 


(3.105) 


Step 5. It remains to estimate the last term of (3.105). According to the fundamental theorem 
of calculus, we have from (3.80) that 


[ sup(r5-l)^< [ [ 1^, 

Jo R+ JO Jn[{s) 

[ 

Jn[ (s) ^ Jn{ 


< e 


< e 


f I €^cioff § 

Jo Jr+ Jo Jr+ ^ 

lo L 


(3.106) 


Plug (3.106) into (3.105) and choose e > 0 suitable small to obtain (3.75). This completes the 
proof of the lemma. □ 
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We obtain the uniform bound of the i/^-norm of {(j), tp, x) uniformly in time t in the next 
lemma. 


Lemma 3.8. Suppose that the conditions listed in Lemma 3.1 hold. Then 

\\{cp,p^,d)mi + 

0<t<T 


sup \\[(j),li),d){t)\\l+ ( + ||(V'a;,l?x)(i)||? 

<t<T Jo 


At < cl 


(3.107) 


where the positive constant C 2 depends only on inf 3 ^gR^{po( 2 ;), ^* 0 ( 2 ^)} and ||(^o, V'o,'do)||i- 
Proof. First, plugging (3.74) into (3.45), we deduce 


/R+ 


^ 1 + 

>0 JR+ Jo JR+ 


'4’x+'&x + 


02 


< 1 , 


(3.108) 


where we employed (3.9) and (3.75) in the last inequality. 
Next, multiply ( 3 . 1)2 by 'ipxx/P to derive 


2 


fjttpx + ^Utpl 


, 1 .. .,.2 , ..^xx {P-P)x_, / 2 _, 

H" n^x'4^x M — '4^xx '4^x 

2 p p p 


Integrating this last identity over (0,t) x K.+ , we obtain from (3.74) and Cauchy’s inequality that 


/ -^lAx- f V’^(s,0)ds+ f [ 1: 
JRo- JO JO JR+ 


2 

XX 


/R+ 

< 1 


(3.109) 


{P-P)l + fl+\Ux\i;l + \f;x 


10 JR+ 

Apply Sobolev’s inequality and (3.75) to obtain 


Pl{s,0)ds+ f [ \Px\^< /‘||V’.||||V’..|I+ f UxW^lPxxt^ 

Jo Jr+ Jo Jo 

f j + f\Uxf + Ux\\- 

Jo JRo- JO 


< 


< 


C(e) 


/o 

1+ sup ||V’x(s)||^ 

0<s<t 


(3.110) 


Pi 


/O JR+ 

Using {P — P)x = RiPPx + P'dx + P^x + dpx), we derive from (3.83), (3.108), (3.75) that 


{P-P)l+\Ux\Pl 


< 


\{0px,px,"dx)p + \{px,Ux,0x)\'^\Al'‘^ 


(3.111) 


< 


1 + II^I|l“([0,T]xR+)- 


According to (3.83) and (3.84), we have from Lemma 2.2, (3.22), (3.8), (3.9) and (3.108) that 
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We plug (3.110)-(3.112) into (3.109) to get 


sup Ux{tW+[ [ V'L ^ 1 + ll^'IU”([ 0 ,T]xR+) + sup 

)<i<T Jo 0<t<T 


Jo iK_|_ 

Then Young’s inequality yields the estimate 

/o Jm+ 

Next, multiply ( 3 . 1)3 by ’dxxjp and integrate the resulting identity over (0,T) x R+ to have 


sup Ux{t)\\^+[ [ ^ 1 + II^IU”([ 0 .T]xR+)- (3.113) 

0 <t<T Jo JRj. 


f / ^idx + K 


f ^ix Cy 

Jm+ P 2 , 


f 

/o M+ 

which combined with (3.74) implies 

rT 


T5o^dx+ / 
1 + Jo 


Cvwdx - p— + ROtjjx - — 

P P J 


dx 


[ dldx+[ [ C<1+/ [ [u^dl + \\i.x\\l^^^l + e^rx + h^] 

Js.+ Jo Js.+ Jo Js.+ 

<1+ r(i+ii^iiiiv.,ii)ii^.f + ruxfuxxw 

Jo Jo 

+ II^IIl“([0,T]xR+) / / + / / /I- 

Jo Jr+ Jo Jr+ 

From (3.113), we have 

/ UxfUxx\\< sup llV’xf / (llV’xf + ||'0x;rf) < 1 + ||6'|lioo([o.T]xR+)- 
Jo 0 <t<T Jo 

In light of (3.108), (3.113) and (3.75), we then obtain 


sup / 'd^.dx -\- 


0<t<T . 


R-f Jo 

Finally, it follows from (3.75) and (3.115) that 
|2 


'dxx ^ 1 + II^IIl°°([o,t]xR+)- 


11^ “ ^IIl°°([ 0 ,T]xR+) ^ sup ||^5(t)llll^?a:(0ll ^ 1 + II^IIl“([ 0 ,T]xR+)- 


0<t<T 


from which we have 

< 1 for all {t,x) £ [0,T] x IR+. 
Combine (3.109), (3.113) and (3.115) to give 


sup f [(j)l+^l + 'dl]dx+ f f [Oct)l + ^l^ + 'dlJ<l, 

0 <t<TjR, Jo JR, 


(3.114) 


(3.115) 


(3.116) 


which together with (3.75) yields (3.107). This completes the proof of the lemma. 


□ 
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3.4 Local lower bound of temperature 


In this subsection, we employ the maximum principle to get the lower bound for the tempera¬ 
ture, which does depend on the time t. 


Lemma 3.9. Suppose that the conditions listed in Lemma 3.1 hold. Then 

infR+ 0{s, •) 


m^f 6>(t, ) - s) -t 1 


for 0 < s < t < T, 


(3.117) 


where the positive constant C 3 depends solely on infa,gR^{po( 2 ;), 6 * 0 ( 2 ;)} and ||(</' 0 ) V’ 07 ^o)l|i- 
Proof. It follows from ( 1 . 1)3 that 6 satisfies 


, /3 ^ ^ r 2 ^ 

Ut uOx ^xx — 

Cy P Cy P ^ p 


> -- 


p2 


R^p 
4:pCyP 4:pCy 




Hence we deduce from (3.74) that 


6t -\- u9x - 9xx + > 0. 

CvP 

Let 0 := 0 - 0 with 9 := c,iniT+e(s!'.){t-s)+i - We observe 


and 


0|x=o.oo > 0, 0|t=, > 0, 


0t + uQx -0XX + 1 ^ 2(0 + ^)0 — 0t + u9x - 9XX + C20^ > 0. 

Cy P CyP 


Applying the weak maximum principle (see [3, Section 7.1]), we have that 0(t,x) > 0 for 0 < s < 
t < T and x G K.+ . This completes the proof of the lemma. □ 


3.5 Proof of Theorem 2 

This subsection is devoted to proving the stability of the superposition of a rarefaction wave 
and a non-degenerate stationary solution, i.e. Theorem 2. To this end, we first give the local 
existence of solutions to the problem (3.1)-(3.2) in the following proposition. It can be proved by 
the standard iteration method (see [ 10 ] for example) and hence we omit the proof for brevity. 

Proposition 3.10 (Local existence). Suppose that the conditions in Theorem 2 hold. Let M, 
Ai and A 2 be some positive constants such that j|((()o,V'O:'*^o)||i < M, 4)o{x) + p{x) > Ai and 
'9o{x) + 9{x) > X 2 for all x € M+. Then there exists a positive constant Tq = To{Xi, X 2 , M), 
depending only on Ai, A 2 and M, such that the problem (3.1)-(3.2) admits a unique solution 
(</.,V',^?) e A(0,ro;iAi,iA2,2M). 

Next we will give the proof of Theorem 2 in six steps by employing the continuation argument. 
Step 1. Let n and Xi {i = 1,2,3) be some positive constants such that l|(())o,^o,^^o)||i < n and 

00 ( 2 ;) > Ai, 00 ( 2 ;) > A 2 , 9{t,x)>X3 foralH,x>0. 

Set Ti = 128A3^C'|, where C 2 is exactly the same constant as in (3.107). Applying Proposition 
3.10, we infer that the problem (3.1)-(3.2) has a unique solution € A(0,ti; ^Ai, AA 2 , 2 n) 

for some positive constant 

ti = min{Ti,To(Ai, A2,n)}. 
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Let 0 < 5 < (5i with 

“ ^^2,2 n) (5i = Co- 

Then we can apply Lemmas 3.5, 3.8 and 3.9 with T = ti to obtain that for each t G [0,ti], the 
local solution constructed above satisfies that 


6{t,x) > 


X 2 

C 3 A 2 T 1 + 1 


=:C4 


for all X S M+, 


(3.118) 


and 


< p{t, x) < Cl for all x G M+, 


rt r 


10 


+ / \\V9(l)^{s)f + \\{^pa:,da:)is)fi 


ds < C^. 


(3.119) 


Step 2. If we take {(j), ip, '&){ti, •) as the initial data, we can apply Proposition 3.10 and extend the 
local solution {(j),ip,i!}) to the time interval [ 0 ,ti + 12 ] with 


^2 — min{ri — fi, To(Cj^ C4,02)}- 


Moreover, for all {t,x) G [ti,ti + 12] x K+, 

p{t,x)>lc^\ e{t,x)>XC4, \\{cP,pj,d){t)\\i<2C2. 

Take 0 < d < min{(5i,52} with 

S(iC'r\iC4,2C2)(52 =eo. 

Then we can employ Lemmas 3.5, 3.8 and 3.9 with T = ti +12 to deduce that the local solution 
{(j),ip,'&) satisfies (3.118) and (3.119) for each t G [0,ti +t 2 ]- 

Step 3. We repeat the argument in Step 2, to extend our solution {(p, ip, i9) to the time interval 
[0,ti +t 2 + ^ 3 ] with 

ts = min{Ti — (ti + ^ 2 ), 'ro(Ci C 4 , C' 2 )}. 

Assume that 0 < <5 < min{ 5 i,( 52 }- Continuing, after finitely many steps we construct the unique 
solution {(p,ip,i}) existing on [0,Ti] and satisfying (3.118) and (3.119) for each t G [0,Ti]. 

Step 4. Since Ti > I 28 A 3 and 


sup \\d{t)\\l+ f \\d^{t)\\ldt<Ci, 

0<t<Ti JTi/2 

we can find a tj, G [ri/2,Ti] such that 

ll^(i'o)ll < C2, \\Mt'o)\\<kc-^xl 

Sobolev’s inequality yields 

||^(i'o)l|L-<^ll^(io)ll^ll^x(t'o)||^<iA3, 

and so 

> 0{tQ,x) - ||'d(to)||L~ > 5 A 3 for all x G R+. 

We note here that 


||(((), V','d)(to)||i < C' 2 , p{t'fj,x) > foralla;GK+. 
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Applying Proposition 3.10 again by taking (cj), ip, •) as the initial data, we see that the problem 
(3.1)-(3.2) admits a unique solution {(p, ip, i9) € A(tg, tg + jAg, 2 C 2 ) with 

t; =min{Ti,To(Cr\iA 3 ,C' 2 )}. 

If we take 0 < i5 < min{di, 52 , dg} with 

S(iCr\iAg,2C'2)5g = eo, 


then we can apply Lemmas 3.5, 3.8 and 3.9 with T = tg + t'l to obtain that for each time t £ 
[to,to + t'l], the local solution {(p,ip,i}) satisfies (3.119) and 


9{t, x) > 


infR^ 6»(t'g,-) 


> 


As 


Ca infR+d(tg, •)’^i + 1 CsAgTi + 2 


=: C 5 for all x £ IR.+ . 


(3.120) 


Step 5. Next if we take {(p,ip,'d){tg + t[,-) as the initial data, we apply Proposition 3.10 and 
construct the solution {cp, ip, d) existing on the time interval [ 0 , tg + tp + ^ 2 ] with 

4 = min{ri - t'l, ro(Cf \ Cg, C 2 )} 


and satisfying 

p(t,x)>icr\ 0{t,x)>^C5, ||(0,V',d)(t)||i <2C2 

for all (t,x) £ [tg + t'ljtg + t'l + ty X K+. Let 0 < 5 < min{5i, 52 ,5g, 54 } with 

S(iCr\iC5, 2 ^ 2)54 = eg. 

Then we can infer from Lemmas 3.5, 3.8 and 3.9 with T = tg+t'i+t^ that the local solution {cp, ip, d) 
satisfies (3.120) and (3.119) for each t £ [tg, tg + t'l + t^j. By assuming 0 < 5 < min{5i, 52,53, 54 }, 
we can repeatedly apply the argument above to extend the local solution to the time interval 
[0,tg + Ti]. Furthermore, we deduce that (3.120) and (3.119) hold for each t £ [tg,tg + Ti]. In view 
of tg + Ti > ITi, we have shown that the problem (3.I)-(3.2) admits a unique solution {(p,ip,d) 
on [0, iTi]. 

Step 6. We take 0 < 5 < min{5i,52,5g, 54 }. As in Steps 4 and 5, we can find tg £ [tg + 
ri/2,tg + Ti] such that the problem (3.I)-(3.2) admits a unique solution {(p,ip,i}) on [0,tg' + Ti], 
which satisfies (3.120) and (3.119) for each t £ [to,to' + Tij. Since tg + Ti > tg + |Ti > 2Ti, 
we have extended the local solution {(p,ip,'d) to [0,2Ti]. Repeating the above procedure, we can 
then extend the solution {<p, ip, d) step by step to a global one provided that 5 < min{5i, 52, 5g, 54 }. 
Choosing e 2 = min{5i, 52,5g, 54 }, we derive that the problem (3.1) has a unique solution {(p, ip, d) £ 
A(0, 00 ; min{C 4 , Cg}, C 2 ) satisfying (3.119) for each t £ [0,oo). 

Therefore, we can find constant Cg depending only on infa;gR^{pg(x), ^ 0 ( 2 ;)} and || [(pg,ipQ, do)|| 1 
such that 

pOO 

sup \\{p,ip,'&){t)\\l+ j [\\p,,{t)\\^+ \\{ip,„d,,){t)\\l]dt<ci, 

0<t<c>o JQ 

from which the large-time behavior (1.24) follows in a standard argument (cf. [23]). This completes 
the proof of Theorem 2. 
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